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■ Abstract 

(N ■ 

We prove that the Schwarzschild black hole is linearly stable under electromagnetic 
, and gravitational perturbations. Our method is to show that for spin s = 1 or s — 2, 

' solutions of the Teukolsky equation with smooth, compactly supported initial data 

, outside the event horizon, decay in L^^. 

^ '■ 1 Introduction 

■ In polar coordinates {t,r,-d,(p), the line element of the Schwarzschild metric is given by 

■ A r2 
O , ds^ = g.k dx^x^ = —dt^ - — dr^ - d^^ - sin^ ^ dif^ 
^•G . r A 

o 



bJO. 



where 



qh. a = - 2Mr 



and M is the mass of the black hole. The zero ri := 2M of A defines the event horizon of 
the black hole. The evolution of a massless wave ^> of general spin s in the Schwarzschild 
geometry is described by the Bardeen-Press equation [5]. Teukolsky [22] later generalized 
^ ' the equation to the Kerr geometry, and therefore the equation is usually referred to as the 

■ Teukolsky equation |6j . We here work with a particularly convenient form of the Teukolsky 

equation due to Whiting [25] : 



drAdr - ^ {r'^dt - (r - M)s)^ - Asrdt 

+5costf sin^ ?? dcos^ + . \ „ {d^ + is cos -df 
sm n 



Ht,r,^,ip) = 0. (1.1) 



This is a second order scalar wave equation, having complex coefficients if s 7^ 0. (Note 
that our wave function differs from the function ip in [22] by a power of A, <I> = A^ip.) 
The parameter s is either integer or half-integer valued. The case s = gives the scalar 
wave equation. Of particular physical interest are the cases s = |, 1, 2, which correspond 
respectively to the massless Dirac equation, Maxwell's equations and the equations for 
linearized gravitational waves. 
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In this paper, we consider the Cauchy problem for the Teukolsky equation with 
initial data 

^\t=o = ^0 , dt^\t=o = ^1 , (1-2) 

which is smooth and compactly supported outside the event horizon. Our main theorem is 
the first rigorous result on time-dependent solutions of the Teukolsky equation for higher 
spin, and proves linear stability of the Schwarzschild black hole under electromagnetic and 
gravitational perturbations. 

Theorem 1.1 For spin s = 1 or s = 2, the solution of the Cauchy problem hl.li 
for (<I>o,<J>i) G C^((ri,oo) x 5^)^ decays in L^^((ri,oo) x 5^) as t ^ -oo. 

The study of linear stability of the Schwarzschild geometry was initiated in 1957 by Regge 
and Wheeler [21|, who discussed mode stability for metric perturbations. In the case s = 0, 
the Cauchy problem (jl.H II. 2p was considered (for more general initial data) by Kay and 
Wald [IS], and they obtained a time independent L°°-bound. Decay in was proved 
in [T2] in the Kerr geometry, and worked out in [16j in the Schwarzschild geometry. 
Related results for s = in the Schwarzschild geometry were obtained in [201 18] . If s = ^ , 
local decay was proved in the Kerr geometry in [9j (for both the massive and massless 
case), and an exact decay rate was given in the massive case [TOj. Up to now, for higher 
spin s = 1 (Maxwell's equations) and s = 2 (linearized gravitational waves) only mode 
analyses have been carried out for the Teukolsky equation; see [19| for a numerical study 
and [25] for a rigorous proof of mode stability. We also mention the two papers on the 
Regge- Wheeler equation [Hill]. In the first paper pointwise decay of solutions is proved, 
whereas in the second paper time integral of solutions are estimated locally in space. 

To consider the limit t — > — oo (and not t — -l-oo) is purely a matter of convenience. 
To see this, first note that in a general space-time, a massless field of spin s ^ satisfies 
a coupled system of 2s + 1 complex, first order partial differential equations. As shown 
by Teukolsky [22], this system can be decoupled in the Kerr geometry by multiplying 
with a suitable first order differential operator. Then the first component of the system 
satisfies the Teukolsky equation (jl.ip . whereas the last component also satisfies p.lj) . 
but with s replaced by —s. From either the first or the last component, all the other 
components can be obtained by applying the so-called Teukolsky-Starobinsky identities, 
see [6]. In view of this, we may restrict attention to the Teukolsky equation (jl.ip for 
either s or —s. Next, we point out that the Teukolsky equation (jl.ip is invariant under 
the transformations {t,s,i!},(p) {—t,—s,'&,—ip). We thus see that Theorem 11.11 also 
makes a similar statement on the solution $ of the Teukolsky equation for spin — s in 
the limit t +oo. Since the Teukolsky-Starobinsky identitities relate the solutions of 
the Teukolsky equations for ±s to each other, obtaining decay for the spin —s equation 
as t — > oo immediately yields decay for the spin s equation as i — > oo. Thus there is no 
loss in generality to consider in Theorem 1 1.1 1 the case i — > — oo. This case will turn out to 
be most convenient if one uses the sign conventions in [22] as well as the factor e~*^* in 
the time separation. 

Let us specify in which sense the Teukolsky equation governs the linear perturbations 
of the Schwarzschild black hole. For spin s = 1, the Teukolsky equation takes into account 
all electromagnetic perturbations except for adding a constant electric charge (see [22\ p. 
644]), thus perturbing Schwarzschild to the Reissner-Nordstrom space-time. For s = 2, the 
Teukolsky equation describes perturbations of the Weyl tensor, and it is a quite difficult 
task to reconstruct metric perturbations from a solution of the Teukolsky equation; for 
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details see [TH', '26] . It is important to keep in mind that the Teukolsky equation excludes 
perturbations of Schwarzschild to the Kerr space-time, but does take into account all 
other regular perturbations (see [22] )• Hence our Theorem shows linear stability of the 
Schwarzschild black hole under all perturbations, except for linear perturbations to the 
stationary Kerr-Newman black hole. 

We now briefly discuss energy conservation and its role in our proof. In the Schwarz- 
schild geometry, the physical energy can be written as the spatial integral of a positive 
energy density. More precisely, in the cases s = and 1, this energy is obtained by inte- 
grating the normal component of the vector field T*q , where Tij is the energy- momentum 
tensor corresponding to the spin s field, 

E = [ T.,,u' = I T%dfi, 

it=const X'^'' / Jt=const 

where i' is the future-directed normal and dfi is the integration measure on the hyper- 
surface t = const. This energy is conserved because the energy- momentum tensor is 
divergence-free and dt is a Killing field. Using the dominant energy condition and the fact 
that dt is timelike, it is easy to verify that the energy density is indeed non-negative. In 
the case s = 2, a conserved energy of the gravitational field is given by the integral of the 
Bel- Robinson tensor Q (see for example ^17j p. 42ff]) 

E = [ Q\oodfi. (1.3) 

This energy is the sum of the gravitational energy of the Schwarzschild metric and the 
energy of the gravitational wave. Thus the energy Egy of the gravitational wave is obtained 
by subtracting the energy of the background, 

i^gv = / iQ\oo-iQ^)\oo)dfi, (1.4) 

where is the Bel-Robinson tensor in Schwarzschild. It is shown in the appendix that 
Egy is conserved and quadratic in the perturbation of the Weyl tensor, and that the 
integrand in p.4|) is non- negative. 

For computing the energy density in (jl.4p one needs to know all the components of the 
spin s field. In order to compute all these component functions from a given solution <I> 
of the Teukolsky equation (jl.ip . one can proceed in two essentially different ways. The 
first method is to take <I> as the first component, and to compute all the other components 
inductively using the first order system of differential equations. The second, probably 
more elegant method is to regard <I> as the so-called Debye potential (a generalization of 
the classical Hertz potential in electrodynamics), from which all the field components are 
obtained by differentiation; this is worked out in detail in [3 [24]. In any case, substituting 
the resulting expressions for the field components into ()1.4p . unfortunately for s = 1 or 2 
one gets very complicated expressions involving higher derivatives of which seem very 
difficult to handle. For this reason, we are unable to use the explicit form of the energy 
density. In particular, we cannot work with a corresponding energy scalar product. As 
a consequence, the associated Hamiltonian will not be a symmetric operator, and thus 
we cannot use spectral theory in Hilbert spaces. The main technical difficulty of the 
present paper is to prove completeness and decay without using the spectral theorem. 
Nevertheless, we will make use of the existence of a positive energy density a few times, 
without refering to its explicit form. 
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The main step in the proof is to derive an integral representation for the propagator, 
whereby we integrate over the real line M together with another line parallel to M (see 
Theorem 18. 5p . The latter line integral reflects the fact that the essential spectrum of 
the Hamiltonian has a contribution in the complex plane. The appearance of a complex 
essential spectrum can be understood from the fact that the Teukolsky equation (jl.ip 
involves first order time derivatives, which after time separation e~*"* with real uj lead 
to complex potentials in the resulting radial equation. These complex potentials make 
it impossible to apply standard techniques used for 1-dimensional Schrodinger equations. 
In particular, the fundamental solutions of the radial ODE behave asymptotically near 
infinity like a power of r times a plane wave (and no longer just like plane waves as 
in the case s = 0). This requires us to develop new techniques like obtaining refined 
WKB estimates for Jost solutions with complex potentials, and working with non-closed 
integration contours. To prove completeness, we use an idea in Bachelot [3], which reduces 
the completeness problem to obtaining resolvent estimates for large values of the spectral 
parameter lo. 

We remark that in the case s = ^ , Theorem 18.51 gives an integral representation for 
the propagator of the massless Dirac equation, which is considerably different from the 
integral representation obtained in [9]. This surprising fact will be discussed in Section [TOl 

This paper is organized as follows. In Section [2] we separate out the time and angular 
dependence in the Teukolsky equation and obtain the radial ODE. In Section [3] we con- 
struct holomorphic families of Jost solutions of the radial equation which have prescribed 
asymptotics near the event horizon or at infinity. In Section U] we write the Teukolsky 
equation in Hamiltonian form and express the resolvent of the Hamiltonian in terms of 
the Jost functions. Section [5] is devoted to the derivation of WKB estimates, which give 
precise bounds for the Jost functions asymptotically near the event horizon and at infinity, 
and also globally if |a;| is sufficiently large. Using these estimates, in Section [6] we study 
the decay properties of the resolvent for large values of the spectral parameter. These 
resolvent estimates allow us in Section [7] to express the propagator in terms of contour 
integrals, thereby also obtaining a completeness result. In Section [8] we use classical Whit- 
taker functions together with an energy argument to show that the integration contour 
can be deformed onto both the real axis and another line parallel to it. In Section [9] we 
prove Theorem 11.11 using a Riemann-Lebesgue argument for a finite number of angular 
modes, together with an estimate for the remaining infinite number of modes. Finally, 
Section [To] is devoted to general remarks on our integral representation in the case s = ^ 
and on possible extensions of our methods to the Kerr geometry. 

2 Separation of Variables 

Using spherical symmetry, we can separate out the angular dependence with the usual 
multiplicative ansatz 

<!>{t,r,'d,ip) = R{t,r) Y{^,ip) . 

The spin-weighted spherical harmonics Y = gYim with / = |s| + 1, . . ., m = —I, —I + 
1, . . . , / (see p[4j) form an eigenvector basis of the angular operator 

A = -dcos-dsm^ dcoss r4-: (^ip + is cos , 

sm V 

on L^(5^), corresponding to the eigenvalues A; = l{l -|- 1) — (note that our angular 
operator is related to the operator in [E] by ^ = — — s^)- Restricting attention to 
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one angular momentum mode, the Teukolsky equation reduces to 



1 

A 



dr/^dr - ^ (r^dt - (r 



M)s) -isrdt - A 



R{t,r) = 



(2.1) 



where we set A = A; . We transform to the Regge- Wheeler variable z/ G M defined by 

— = — , so M = r + 2M logfr - 2M) , 
dr A 

which maps the event horizon r = 2M to u = —oo. Furthermore, setting 

(j){t, r) = r R{t, r) , 

the Teukolsky equation becomes 



(2.2) 



^dur^du- - 4 (r^dt - (r - M)sf - Asrdt - A 
A r A 



(t>{t,r) 



Applying the identity duv'^du = rd'^r — r{d'^r), we can write this equation in the simpler 
form 



dl - a 



(r - M)s 



4sA dlr A 
— 3- Ot A — 



<P{t,r) 



0. 



(2.3) 



Using the time translation symmetry, we can further separate out the time dependence 
with the ansatz 

0(t,r) 



Then the Teukolsky equation reduces to the ODE in Schrodinger form 

- ^ Hu) + Viu) <P{u) = , 



where the potential V is given by 

"2(r-M) 4A 



Viu) 



-co + ISUl 



ir-Mfs'^ dlr , A 
+ - + ^ + A^ 



(2.4) 



(2.5) 



(2.6) 



Note that in the case s ^ V is complex even for real iv. 



3 Construction of the Jost Solutions 

In this section we construct Jost solutions and (j) of the Schrodinger equation (j2.5p . 
which are defined by their asymptotic behavior near the event horizon and near infinity, 
respectively. Near the event horizon, the potential has the limit 



lim V{u) = — where := 

Writing the Schrodinger equation in the form 

{-dl -n^)<t> = -w{u) <t>iu) , 



IS 

4M 
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the potential W behaves near the event horizon hnearly in (r— ri), and thus has exponential 
decay in the Regge- Wheeler coordinate ()2.2p near n = — oo. More precisely, for u near — oo 
there is a constant c > such that 

|VF(n)| < c e2i? . 

Using this exponential decay, (p can be constructed exactly as in [121 Theorem 3.1] (see 
also [21 IS]). The properties of (p are summarized in the following theorem. 



Theorem 3.1 For every iv in the domain 

D_ = 



UJ 



there is a solution (/>_ of 112.5]) having the asymptotics 

-iQu 1 /,\ 



lim e 



1 



It— > — OO 



lim (e-^^"(^_(u) 



. 



(3.1) 



These solutions can be chosen to form a holomorphic family, in the sense that for every u G 
M, the function (u) is holomorphic in u) ^ D_ . Similarly, on the domain 



ijj 



1 



A 



4M 4MJ 

there is a holomorphic family of solutions 4>+ of 112. 5\) with the asymptotics 



lim e 

u—>—oo 



<j>+{u) 



lim ( e 

u— >— oo 



iQu 1 



<j)+{u) 







Near infinity, the potential has the following asymptotic form, 



N 9 2isio /logn , 

V{u) = -uj^ - — + O { -|- ) . (3.2) 



In the remainder of this section we always assume that u ^ 1. In the case s = 0, 
the solutions cj) were constructed in [121 IE] ! thus we assume in what follows that s > ^ . 
Because of the non-integrable M~^-term in V, the standard Jost solution method [2] cannot 
be implemented. We choose uq so large that V has no zeros on [uq, oo). We introduce the 
WKB wave functions d and d by 

a(n) = cy(n)-3exp ( / VV\ , d{u) = cV{u)-^ eicp (- f Vv\ , (3.3) 



with constants c, c 7^ to be determined later. To explain the sign convention for yV, we 
first note that taking the square root of (13. 2p gives 



V^=±(..-l)+0(l5^). (3.4) 



Our sign convention is 

+ if Imw < 



.„ ^ (3.5) 

iflma;>0. ^ ' 
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Thus if Imw / 0, the real part of ^yV{u) is positive for large u and so a decays at plus 
infinity. Furthermore, we note that our sign convention does not change if to approaches 
the real line from below. Also, we point out that for real u, the function a does not decay 
at infinity, but increases polynomially like n*. 

The functions d and d are solutions of the equation 



La = , 

where L is the differential operator defined by 



(3.6) 



L 



-dl + Vo and Vq := V 



V" 



+ 



4V16\V 



V 



Writing the Schrodinger equation (|2.5p as 

= -W(p, 

we see that the new potential W := V — Vq is integrable, since 



\Wiu)\ < 



1 



1 + \uj\ u 



on [uq, oo) 



(3.7) 



(3. 



Since the WKB wave functions d and d form a fundamental system for ()3.6p . we can 
use them to construct a Green's function for the operator L. In what follows Q denotes 
the usual Heaviside function. 



Lemma 3.2 Under the sign convention \3. 5|) . the function 



S{u, v) 



Q{v-u) {V{u)V{v)y 



exp 



exp 



is, for all u,v > uq, a distributional solution of the equation 

LuS{u,v) = 5{u — v) . 

Proof. We make the ansatz 

S{u,v) = @{v — u) {ci{v) d{u) + C2{v) d{u)) 
and determine the coefficients ci and C2 from the conditions 

0, lim duS{u,v) = —1. 



lim S{u, v) 

u /'v 



This gives (13. 9p . and a straightforward calculation yields (I3.10p . 



(3.9) 



(3.10) 



We now make the perturbation ansatz 

oo 



n=0 



where the 0^"^ are defined by the iteration scheme 

,/,(0)(n) = d(n) 

/>oo 

(n) = -/ S{u,v)W{v)<j)^^\v)dv . 



(3.11) 



(3.12) 
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Before stating the next theorem, we must study the asymptotics of a near infinity. Car- 
rying out the integral in (j3.3p using (j3.4p , we obtain 

a(n) ~ gWi^iog" = ^±s^Ti^u^ 
Due to our sign convention (|3.5p , we find that 



~ I n^e— iflma;<0. ^^"^^^ 

We next prove that the perturbation series (13. lip converges to a solution cp of the full 
equation (|3.7p having the same asymptotics as a. 

Theorem 3.3 On the domain := {tu \ tu ^ and Imcj > 0}, there is a family of 
solutions of ([g. 7[j, holomorphic in the interior of E^, having the asymptotics 



lim e"*^" = 1 , lim e"*^" </>+(«)) = . (3.14) 



Likewise, on the domain E- := {u | tj 7^ and Imu < 0}, there is a family of solu- 
tions (p- (u) of Ii3. 7\ ), holomorphic in the interior of E^ , with the asymptotics 



1, lim fn-"e*'^"0_(u)y = 0. (3.15) 

M— +00 V / 

Proof. From the definitions (j3.31 13. 9p and our sign convention (]3.5p , it is obvious that 
\a{u)\ < d\V{u)\-i exp(^- j \ReVv\^ 



1 



\S{u,v)\ < Q{v-u)\V{u)V{v)\-i |ReVF| 
where we take 

d = exj)(^2 j ' \ReVv\^ , 

with ui chosen so large that the real part of the square root of V is positive for all u > ui. 
We will show inductively that for u > uq, 

(t>^^Hu)\ < cd^ \V{u)\~'^ e^p(^- J iRe^l^ , (3.16) 

where 

C = " 

2(l + |w|) 

and c is as in ()3.8p . The case / = is obvious. Assume that (j3.16p holds for a given /. 
Then 



|0('+^)(^)| < hd^\V{u)\-\e-^^-o\'^^\ r -^^^dv 

L. i -|- 10^1 V 

= ^' n^('+i)(^ + l)! '^^"^' ' ■ 
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The estimate (|3.16p shows that the series (|3.1ip converges absolutely, uniformly for u > 
uq. Similarly, one can show that the series obtained by diffentiating (j3.1ip termwise again 
converges in the same sense. Thus we can differentiate the series termwise, thereby showing 
that (/> is a solution of (j3.7p . According to (j3.13p , we can choose c such that the function (j)^^'^ 
satisfies the boundary conditions (j3.14p or (j3.15p . and since the estimate (j3.16p involves 
a factor u~'^\ it is obvious that (j) also satisfies the first relation in (|3.14p or (|3.15p . The 
second relations are obtained by differentiating (|3.1ip and (j3.12p with respect to u; a 
lengthy but straightforward calculation yields the result. 

To prove analyticity in lo, we first note that a, W and 5'(n, v) are obviously analytic. 
Hence each c/i^'^ is analytic being an integral of analytic functions. Since the constants c 
and d, C in (j3.16p can be chosen locally uniformly in a;, we conclude from Morera's theo- 
rem that (j) is also analytic. ■ 



4 Hamiltonian Formulation, Construction of the Resolvent 



At this stage, we do not know whether the separation ansatz (j2.4p will give us a complete 
set of solutions of the time-dependent equation (j2.3p . To remedy this situation, we shall 
write the equation in Hamiltonian form. To this end we set 



idt(j){t,r] 



and obtain 



and 



with 



H 



1 

a [5 



a 



P 



2 , (r-M)2s2 



-d,+ 



A 



IS 



2(r - M) 



4 A 



(4.1) 

(4.2) 
(4.3) 



The Hamiltonian H can be considered as an operator on the Hilbert space 



n 



H^'\R,du)eL\R,du) , 



densely defined on the domain T>{H) = 5(M)^, the Schwartz functions. Note that H is 
not symmetric on Ti. 

We assume for the rest of this section that 



Imw 



0, 



4M 



(4.4) 



For a given u; satisfying these conditions, we will show that the resolvent of H exists, and 
we will express it in terms of the Jost solutions. Depending on the sign of Imw, we let 
be the function </>+ or (j)- of Theorem 13.31 respectively. If a; G I?+ n D- , there are two 
Jost solutions (j)± near the event horizon, one of which decays exponentially, the other of 
which grows exponentially. We choose the solution with exponential decay. Thus in the 
case Imu; > s/(2M) we let (f) = (/>+, whereas in the case Ima; < s/(2M) we let <p = <j)-. 
The next lemma relies crucially on Whiting's mode stability result |25j . 
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Lemma 4.1 For any lo satisfying \4-4% the Wronskian 



w{(p,(p) := (P'^-(P^' (4.5) 



IS non-zero. 



Proof. If the Wronskian were zero, the solutions cf) and (j) would be linearly dependent. 
Then there would be a solution (p decaying exponentially fast at both u = ±00. If Imo; < 0, 
such solutions have been ruled out by Whiting (Note that Whiting considers the 

case s < in the limit t — > +00. Using the symmetries {t, s, ip) — > (— t, —s, —^p), this 
is equivalent to considering the case s > in the limit t — > —00, which corresponds to 
mode solutions in the lower half plane as considered here). 

In the case Im uj > s/(4M) and Re a; 7^ 0, we use that (/> is a solution of ()2.5p to obtain 

Using the exponential decay of </) as u ^ ±00, we can integrate by parts to get 

= -2 (0,(Imy)<A)^.(j,). (4.6) 
On the other hand, we see from ()2.6p that 

"2(r - M) 4A" 



s 

ImV = — 2Reu; Imw — — 

2 



(4.7) 



A short computation shows that the round bracket is strictly positive. Thus Im V is either 
strictly positive or strictly negative, contradicting ()4.6p . 

In the final case Imw > s/(4M) and Reu; = 0, we see from ()4.7p that V is real, and a 
short computation using (12. 6p shows that it is even strictly positive. Using that according 
to (j3.ip . the fundamental solution (/> is positive and increasing near u = —00, we conclude 
that (j) is convex. Hence it cannot be a multiple of the function (j), which decays at infinity 
according to (I3.14p . ■ 

This lemma allows us to introduce the Green's function G{u, v) of the Schrodinger equa- 
tion (|2.5p by the standard formula 

G{u,v) = { ^) i ') i ~ ■ (4-8) 

■w{(j),(p) [ 4>{u)(p{v) itv<u. 

It satisfies the distributional equation 

^ + V{u)^G{u,v) = 6{u,v). (4.9) 

We let G denote the corresponding operator with integral kernel G{u,v). 

Lemma 4.2 For every uj satisfying G is a hounded linear operator from 

to H^'^{R), and maps C^(M) to S{R). 
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Proof. We restrict attention to the case Imw < 0, since the other case is analogous. To 
prove the first part, we let ip be in L^(R). Then the function Gip can be written as 

{G'4j){u) = ^^[4>(u) / !i){v)^{v)dv + (j){u) / (t>{v) ip{v) dv 



w{(j),4>) 

We consider only the first term, because the second term can be treated similarly. Thus 
our task is to bound the function 

/•oo 

f{u) := <j){u) / ^^{v)^{v)dv (4.10) 



in H^'^. From Theorem l3.1l we know that the solution behaves near the event horizon like 
(j) ~ e*^". Integrating the Wronskian equation (j)'h — 4>h' = 1 via the method of variation 
of constants, we obtain another fundamental solution 

/■o 1 

h{u) = (l){u) / -j—dx. 

Using the asymptotics of 0, one sees that h is bounded near the event horizon by a 
multiple of |e~* |. We conclude that the function (f), being a linear combination of these 
two fundamental solutions, satisfies the inequality 

Using similar arguments near infinity, we conclude from Theorem 13.31 that 

|0(n)| < C7^i^el^°''^l" foru»0. 

Combining these inequalities with Theorems 13.11 and 13. 3^ we have estimates for both (p 
and (j) at both asymptotic ends. Since on any compact set, the solutions can be bounded 
using simple Gronwall estimates, one sees that choosing e = min(|Ima;|, |ImO|), we have 
the following estimate for sufficiently large c, 

ku) Hv) < c e"^ for all v > u . (4.11) 

This estimate gives the pointwise bound 

/•oo 

\f{u)\ < c e-^(^-") |V(t^)| dt;. (4.12) 

J u 

Setting 

g{x) = c6(x) e"*^^ , 

we can write the right of side of (j4.12p as the convolution 5* | V"! • Then using the Plancherel 
Theorem together with the fact that convolution in position space corresponds to multi- 
plication in momentum space, we have 



I2 < ||5'*IV'll|2 = ll^-IV-llb < II5II00 ||''/^||2 

The function g can be computed (ignoring factors of 27r) to be 



00 



g{k) = c I e e*''^ dx - 



e — ik ^ 
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and thus ^ is a bounded function. We conclude that there is a constant C such that 



\2 < C \HfJ\\2 

To get a similar L^-bound on /', we first differentiate (j4.10p . 

POO 

f'{u) := —cl){u)(j){u)ijj{u)+^'{u) I (j){v) ■iIj{v) dv 



Using (j4.1ip . we see that the first term is in L^. To bound the second term, we solve the 
Wronskian equation (14. 5p for (/)' and use the above inequalities to obtain 

|0'(n)| < C^/^el^'^'^l" foru»0. 

In view of Theorem [3TT1 we have similar inequalities for (p' as for (j), and thus we can repeat 
the above arguments with <j) replaced by <p' to obtain 

Wf'h < cuh. 

We conclude that G is a bounded operator from to H^''^. 

It remains to show that G maps Cq^ into the Schwartz class. By iteratively taking the 
derivatives (du + d^) of (j4.8p . we see that {du + dv)'^G{u,v) is continuous in both variables. 
Since for any iIj G Cq°, 

poo POO roo 

du G{u,v)^{v)dv = / {{du + dMu,v))il>{v)dv + / G{u,v)i^'{v)dv , 



where the last integral was obtained by partial integration, it follows that is in C^. 
The higher regularity follows by induction. To prove that G'ljj has rapid decay, we choose u 
to the left of the support of ■;/'• Then 

(G^)(n) = m r^Mdv. 

J-oo w{4),(t)) 

Since ^ has exponential decay, it follows that GV' has rapid decay at u = -co. A similar 
argument using (j) shows that Gip has rapid decay at +00. Differentiating through the 
Schrodinger equation (j2.5p . one sees that the derivatives of 4> and (p also have rapid decay 
at their respective asymptotic ends, implying that all derivatives of Gifj have rapid decay. ■ 



We now express the resolvent of H in terms of G. 

Theorem 4.3 Every complex number to satisfying 1^4 ■4\ ) ^^6-5 the resolvent set of the 
operator H. The resolvent := {H — oj)^"^ has the integral kernel representation 

/oo 
RUu,v)^iv)dv , (4.13) 
-00 

where 

RUu,v) = ')+Giu,v)( M , (4.14) 

\ d{u, V) J \uj {(jj — p(v)) uj J 

and P is defined as in 
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Proof. A short calculation using (j4.1l 14.141 14. 9|) shows that 



{H - uj) R^{u,v) = \5{u-v). (4.15) 

Using Lemma|121 we can use ()4.13p to define as a bounded operator from Ti = i7^'^©L^ 
to itself. 

Let us show that the image of the operator {H — lo) (with domain of definition 'D{H) = 
5(M)^)) is dense in TL. To this end, for given ^ G we choose a sequence £ Cq^ 
with ^ va. 7i. According to Lemma l4.2|, the functions := Ru)^n are Schwartz 

functions. Hence the are in the domain of and from (|4.15p we see that {H — to)^^ = 

We conclude that oj lies in the resolvent set of H and that [H — uj)^^ = R^. ■ 



We end this section by showing that the boundary of the set (14. 4p lies in the essential 
spectrum of H. 



Proposition 4.4 



aessiH) D M U ( M + ^ 



Proof. Let u; G M U (M + j^) and set k = Reuj. We choose a positive test function rj G 
Cq°((— 2,2)) with = 1 and consider for any L 7^ the "wave packet" 

of momentum k, localized in the interval [L^ — 2L^, + 2L^]. A scaling argument shows 
that ||^'k,w,l||l2 = lkllL2) and thus the Hilbert space norm is bounded away 

from zero as L ^ ±00. Furthermore, moving the wave packet to infinity and to the event 
horizon, respectively, we can use the asymptotic form of the Hamiltonian to obtain 

lim \\{H -uj)^^,^^L\\n = ifu; = K 

L— *oo 

ts 

lim WiH -u)^^^^^l\\h = ifw = K + — . 
L^— 00 4M 

Hence uj lies in the approximate point spectrum. ■ 



5 WKB Estimates 

In this section we again assume that ()4.4p is satisfied and that for a suitable constant K > 1 
(to be determined later) one of the following two conditions holds: 

(CI) |u;| >KandueR. 

K 

(C2) a; / and \u\ > —r. 
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By choosing K sufficiently large, we can clearly arrange that the potential V in (|2.6p 
has no zeros. Then the WKB functions d and a are defined by (|3.3p . We choose the 
normalization constants c, c such that 



lim e^""d 

s —iu)u 



lim u e 



a 



lim e-*^"d 



hm u~' e*^" d 



if Im u) > 



4M' 



if Imw < 0. 



The next theorem shows that for large the fundamental solutions <j> and (p constructed 
in the previous section are well- approximated by the WKB solutions. We restrict attention 
to the physically interesting cases s = ^,1,2 (although the method works for arbitrary s 
just as well). We let p be the function 



p{u) = ^/l + u"^ 



and introduce the constant p by 



1 in case (CI) 
p{K/\uj\) in case (C2) . 



Note that in both cases, p{u) > p and \uj\p > K. 

Theorem 5.1 Let s S {^,1,2}. Then there are constants C,K > such that for all lo 
and u satisfying and either (CI) or (C2), the following inequalities hold 



a 



+ 



a 





AC 










1 


< 


and 




+ 






\uj\ p 




a 




a 



< 



4C 
|u;| p 



Proof. We only give the estimate for (j), because 4> can be treated similarly after replacing u 
by —u. By choosing K sufficiently large, we can arrange that for n = 1, 2, 3, 



\Wiu)\ < 



\V{u)\ > 



UJ 



4 ' ■ ■ p 

Furthermore, it follows from our sign convention (j3.5p that for all lo satisfying (|4.4p 



(5.1) 
(5.2) 



ReW > -- + 0(p-^) 
P 



and integrating gives the bound 



■f:2Vv 



< c 



p{u) 



2s 



for all X > u. 



Using ([331EJ1E12]), the functions defined by 



^«(u) 
a{u) 



(5.3) 



(5.4) 
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satisfy the relations 



^(0) = 1 



W{x) 



|l_e-2rVy|^(0(^)rf^. 



(5.5) 



We begin with the case s = 1. We shall prove inductively that for sufficiently large C 
there are constants a^'^ and 6*^') such that for all / > the following inequalities hold, 



u) — a 



(0 



< 



p{u) 



with 



,(0| 



C 
\uj\ p 



(5.6) 



To satisfy these conditions in the case Z = 0, we simply set a^^^ = 1 and b^^^ = 0. Thus 
assume that ()5.6p holds for a given I. Since p >1, ()5.6p implies that 



(5.7) 



The estimates 



W{x) 



< 



< 



c 



dx 



cC^ 



(I^Ip) 



i+i 



1 



— oo 1 + 3; 



dx 



CTT C' 



(klp) 



l+l 



give us control of the first term in the curly brackets in (jS.Sp . To estimate the second term 
in the curly brackets, we first consider the error term in (|5.6p . 



u 2^/V{x) 



p{x) 



< 



< 



1 



\^\ju p{xf \p{u)J p{x) {\uj\py 

C^cc 



dx 



(|a;|p)'+i7n P{xy 
For the constant term in (j5.6p we can integrate by parts 



°° 1 , C^TTCC 

dx < 



(Hp) 



i+i 



u 2,/V{x) 



,(0 



W{x) d ^^,2CVy 
AV{x) dx 



dx 



to get 



4y(u) 



« 2yF(^ 



4y(a;) J 



(5.8) 



< 



c C 
+ 



7 /'OO 



CC 



p{x) 

2|w|'+2p^ p(-u)3 ' (|a;|p)' y„ |w|2p(x)4 V/o(n) 



+ 



C'ccvr 



< 



+ 



C^CC-K 



2|t^|'+2p^ p(-u)3 |wp+2p'p(u)2 - 2(|cj|p)'+l (|c^|p)'+l 
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Choosing (7 > ctt + ccvr + (c/2 + vrcc), the induction step is thereby complete, so that (|5.6p 
holds. 

We choose K > AC. Then the inequality \uj\p > K implies that K/{\uj\p) < and 
thus 



f -1 



15M 



53} - E3 - / 

1=1 1=1 



< 



E 



|a;| p 



C 



< 



AC 



\io\p — C 3|u;|p 



To treat the derivative of 0, we first note that 

(<^(0)' 



a 



Differentiating (jS.Sp and using (j3.3p . we find that 







^1 /"OO 










AV2. 


J u 



W{x) 

In x/yR 



e-2/:v^i?(0(^)Ja; 



Using (j5.2p . one sees that in both cases (CI) and (C2), the square bracket is uniformly 
bounded away from zero, and the integral can be estimated exactly as before. This shows 
that 

OO ^ 

a 



E 

1=1 



{i+i)y !± 



can be made arbitrarily small by choosing K large enough. Thus 

AC 



1 



a' 



< 



?>\u}\p 



and this concludes the proof of the theorem in the case s = 1. 

In the case s = ^ the proof is even easier since we do not need to integrate by parts 
in ()5.8p . Finally, if s = 2, we again consider solutions of (j5.5p . but we replace (|5.6p by 



^«(m) - a 



(0 



-.(0 



6(0 



p{u) p{uf 



< 



p{uf 



(5.9) 



and our task is to prove inductively that there are constants C and K such that for all I 
and uj satisfying the conditions (CI) or (C2), 



,(0 



c 

uj\ p 



(5.10) 



Once these inequalities are proved, the theorem follows as in the case s = 1. Again 
for I = 0, setting a^^^ = 1 and 6^^^ = c^^^ = = 0, there is nothing to prove. The 
induction step follows as in the case s = 1, however we here need to integrate by parts up 
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to three times. We shall not give all the details, but merely consider the term involving b^^\ 
which is representative of all other terms. After integrating by parts twice, we get 




Carrying out the differentiations, we can take absolute values and estimate term by term 
using (15.11 15. 2p . Possibly after increasing K, we get the desired result. ■ 



6 Resolvent Estimates for Large \lu\ 

In this section we assume again that uj is in the range (j4.4p and that condition (CI) holds, 
so that the WKB estimates of Theorem 15.11 are valid. Our goal is to prove the following 
estimate of the resolvent for large which will play a crucial role in the completeness 
proof. 

Theorem 6.1 For every G C^(M)^ and w £ R there are constants K,C = C{u) > 
such that for all to satisfying ^4-4\ )j 

\UJ\ 

Proof Noting that for any -i/' G Co°(M), 

/oo 
Giu,v) {-d^ + V)^{v) , 
-oo 

a short calculation using ()4.131 l^?T^ allows us to write the resolvent for any ^' G Cq°(M)^ 
as 

= j G{u,v)N^ where N = (^"^ ~^ 
Thus, since is linear in lo, the result will hold if we show that for every ^ £ C(J°(M), 

|(G^)(n)| < (6.1) 
I'^l 

Before giving the proof of (j6.ip . we collect a few properties of the potential V for 
large It is obvious from ()2.6p that there is a constant m > such that for all uj 
satisfying (j4.4p and li e M, 

\V{u)\ > \V'{u)\ < m\u\. (6.2) 
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Furthermore, writing the potential in the form 

V = -uj'^ + LoP + f 
(where /3 and / are independent of to), its square root can be written as 



u 



(5{u) f{u) 



Using our sign convention p.Sp together with the fact that (3 and / are bounded functions, 
we conclude that there is a constant m > such that for all oj satisfying (j4.4|) and u G M, 



Re > 



-m . 



Using 



we have 
1 



u) I (f){v) 'il){v) dv + (f){u) I (j){v) ip{v) dv 



w{(t),(t>) 

We first estimate the Wronskian. From (13.31). we have 



cc 



aa 



a' 



V 

XL 



+ Vv] a 



a 



so 



— — a a — — — aa 

a a a a 



cc — - 



V 

3 

4V2 



+ 1 



cc — 



a' a \ / a a 

Applying Theorem 15. H we obtain, possibly after increasing C, 



2cc 



1 



< 



C 



\UJ\ 



Next we consider the second term in the brackets in ()6.4p . We have 



4>{v) ^(v) dv 



6' dv 



J u (p Ju 



A/2 



1 UV'. 



Squaring the identity 



cj)' 1 I6.611 (/^ 



' a 



and using the equation (j)" = V(j3, we find that 

a' 



hl2 



A/2 



a 



V 



,2 \-i 



+ 



161/3 ) 



(6.3) 
(6.4) 

(6.5) 
(6.6) 

(6.7) 
(6.8) 

(6.9) 



dv 



(6.10) 
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Hence 



A/2 



a' 



a 



+ 



v 



12 



2yl 161^3 



t + t.:tt^I , (6.11) 



2yf 



16F3 ) 



and using Theorem 15.11 together with (j6.2p . we see that (j6.1ip is of order l/|w|. Now we 
can estimate ()6.10p termwise to obtain 



(t){v) tp{v) dv 



< -f^ snp\^\ , 
I'^l JC 



where fC denotes the support of ip. 

The first term in the brackets in (16. 4p can be estimated in a similar way. We thus 
obtain 



|(GV)(n)| < ^ sup- 



1 



I'^l veK. \w{(t>, (p) 

Now from (j3.3p and (|6.3p . we get in the case v > u, 



(v) (p{u)\ Q{u -v) + \(j){u) (l){v)\ Q{v - u) 



{u)4>{v)\ = \cc\\V{u)V{v)\~h-^>''^ < |cc| |y('u)F(w)rie™(^-^) < 



cm 



CO 



where C clearly depends on /C. The case u > v can be treated in a similar way. We finally 
apply ()6.9p and possibly increase K to obtain (j6.ip . ■ 



7 An Integral Representation of the Propagator 

In this section we shall express a given G Cq°(M)^ in terms of a contour integral of the 
resolvent. Our method avoids spectral theory and Hilbert space techniques. Instead, it 
uses an idea which we learned from A. Bachelot [3, Proof of Theorem 2.12] and is based 
upon the resolvent estimates of Theorem 16.11 The result in this section is in preparation 
for the integral representation of the propagator which will be derived in Section O 

For given i? > we consider the two contours Ci and C2 in the complex cj-plane 
defined by 

Ci = dBnio) n {Imcu < 0} , C2 = dBnio) n {imcu > ^} , 

both taken with positive orientation; see Figure [H We set C/j = Ci U C2. We can now 
state the following completeness result, valid for any spin s G {^, 1, 2}. 

Theorem 7.1 For every ^' G Cq°(M)^ and n G M, we have the representation 

= -1^ lim / iR^^)iu)doj. (7.1) 

Proof. Since the length of the contour S*! U 5*2 := dBji{0) \ C stays bounded for large R 
(see Figure [T|) , 
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As a consequence, 



lim / — = 1 . (7.2) 

Since our contours (omitting the end points) lie in the resolvent set of H (see Theo- 
rem we know that for every to £ Cr, 

^ = R^{H-uj)-^. 

Dividing by oj and integrating over Cr, we can apply (j7.2p to obtain 



- lim / ^ (R^ (H - io)^){u) . 



'Cr 

But the second term in the curly brackets vanishes in the limit, because using Theorem 16.1 
and the fact that G Co°i we have 



[R^H^>){u) — 



< 



Cr 



C \duj\ ^ 2ttC 
\u}\ ~ R 



Thus (HH) holds. 



Our next objective is to derive an integral representation for the solution of the Cauchy 
problem. We consider the solution ^'(t, u) = {<j),idt(j)) of the separated Teukolsky equa- 
tion ([2T3j) for initial data G C^(M)2. The difficulty is that {R^'^o)iu) only decays in uj 
like 1/io and thus we cannot take the limit i? — > oo in ()7.ip using the Lebesgue dominated 
convergence theorem, nor can we commute differentiation with the limit. To remedy this, 
we derive a finite Laurent expansion of (i?^^'o)(u). 

Lemma 7.2 For every n E N and uo in the resolvent set of H, 

{Ru,^o){u) = 2 h — [RcoH Wo}{u) . (7.3) 
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Figure 2: The contour G^/j = £i U £2- 



In particular, for n = 3 we have for large \uj\ the expansion 



uj + i {uj + iY 

-^o{u) + 2i{H^o){u) + {H^^o){u) 



+ 0{--j . (7.4) 



(uj + i)^ \\u;\^ 
Proof. Dividing the equation Ru){H — u) = 1 by u; gives 

R^-^o = -— + -Ru.H-^o, 

CO CO 

and since H^q is again in , we can iterate this formula to get ()7.3p . The equation (17. 4p 
follows from (j7.3p using the Taylor expansions 

1 1 1 1 i 1 .^/, , 4x 

V oj + i) 



For the application to the Cauchy problem, it is more convenient to deform the con- 
tour Cr to the contour (t/j = U as shown in Figure [2j This contour deformation is 
immediately justified from the analyticity of the resolvent in the respective regions. 

Theorem 7.3 For any spin s G {^,1,2}, the solution of the Cauchy problem for the 
separated Teukolsky equation i2.3\) with initial data "fo = {(p, dtcl))\t=o £ C^(M)^ has the 
following representation: 

= — — lim [ e-''^^ {R^'^o)iu) duj . (7.5) 
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Proof. Let us first verify that the hmit i? ^ oo in (j7.5p exists. To this end, we first 
note that inserting the "counter terms" Ck/{uJ + i)^ does not change the integral in the 
hmit i? — > oo, 

hm / e-"^' {R^^o){u) (ko = lim / e"^'^* ( _ V -^^^ ) . (7.6) 



This is easily verified using the Cauchy integral formula for the closed contour ^rU SiU S2 
together with the fact that the integral over the contours Si and ^2 vanishes as -R ^ 00 due 
to the 0(|a;|~"^)-decay of the counter terms. By choosing the coefficients Ck as in (j7.4p . we 
can arrange that the integrand decays like Thus we can apply Lebesgue's dominated 

convergence theorem to see that the limit i? — > 00 exists. 

Setting t = 0, it follows from Theorem 1 7 . 1 1 that ^ satisfies the correct initial conditions. 
To see that is a solution of the Teukolsky equation, we apply the operator {idt — H) 
to (j7.6p . Since taking the time derivative of the integrand gives a factor of —iuj, whereas 
the WKB-estimates of Theorem 15.11 show that the spatial derivatives of (i?^,^'o)(tt) scale 
like powers of uj, we see that the partial derivatives of the integrand on the right side 
of (|7.6p can all be dominated by a constant times jw]"^. Hence we can interchange the 
differentiations with the limit and the integration to obtain 



(idt-H) lim / e-''^AR^'i/o){u) duj 
= - hm / e-^-^H -co) ({R^^o){u) - Y,-^^] dw . 

If the operator H acts on the factors Ck{u), the resulting expressions are exactly of the 
form as considered after (j7.6p and vanish in the limit. If a; multiplies C2 or C3, the resulting 
terms again vanish in the limit. Hence, using that ci = —^(){u), we obtain 

(idt-H) lim / e-*"^* (i?^^'o)(n) 

= - lim / e-- (Mn)-u;p^)d. = - lim / e^* (P^)du. = 0, 

where we again used the argument after (j7.6p . ■ 



8 Contour Deformations Onto The Real Line 

The objective of this section is to move that part of the contour, which in Theorem 17.11 
lies in the lower half plane, onto the real line. Since by Theorem 13.11 <p is holomorphic 
in a neighborhood of the real line, our task is to analyze for any given G K the Jost 
solutions <j) for uj in the the set 

Ce(a'o) = {l"^ — "^ol < £ and Imw < 0} 

and consider their limiting behavior as w — > wq. We distinguish the cases loq = and ujq 7^ 
0. 
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Ill 

We begin with the first case loq = 0. Qualitatively, if n ^ I^^F i th^ solution (p is well- 
approximated by the WKB solution (see Theorem 15. ip . If on the other hand u <C 
the solution should be close to the solution of the Schrodinger equation for a; = 0. In order 
to match the asymptotics through the intermediate region, we need a separate argument 
based on classical Whittaker functions. 

Lemma 8.1 Let (j) he the J ost functions constructed in Theorem \3. 51 for uj G Ce(0). After 
suitable rescaling, the limit w ^ exists, 

lim Lo'+"^(P = ^0, (8.1) 



where 

a = ^ (^l + 4s2 + 4A-l) . 



The limit function 4>o is a solution of the Schrodinger equation \2. 5|) for a; = and has 
the asymptotics 

r— »oo \ / (2«j I [a + I — s) 

Proof. In order to avoid the difficulties associated with the term proportional to u^"^ \ogu 
in the potential ()3.2p . it is easier to work in the coordinate r. Introducing the function 



^(r) = VAR{r) = — (j){r) , (8.4) 



we can write the Schrodinger equation (|2.5p as 



-^V'(r)+V(r)V(r) = 0, (8.5) 



where the new potential V has the following asymptotics near infinity: 

, 2 isuj + Muj^ s^ + X-2iMsuj - 12]VP uj^ , , 

V(r) = -u?-2 + 5 + 0[r-^). (8.6) 



r 



We first consider the equation (jS.Sp where we simply drop the error term in (18. 6p . Then 
this modified equation can be solved exactly using Whittaker functions [U Chapter 13, pp 
505-508]. The two fundamental solutions are M^^^{z) and M^K,/i(-2), where the parameters 
are given by 

K = s- 2iujM , A* = ^ a/i + 4s2 + 4A - SiMsLO - 48 M^oj'^ , z = 2iujr . 

The function (p clearly is a linear combination of M^^^{z) and Wk_^(2:). Comparing the 
asymptotics for large \z\ [H (13.5.1) and (13.5.2)] with the asymptotics oi (j) (I3.15P . we can 
determine the coefficients of this linear combination to obtain for the function -ip = ^/~K(f)/r 

= (2iw)-^+2^^'^- W.,^{z) . 

Using the asymptotics for small z [H (13.5.6)], we find that again after dropping the error 
term in (j8.6p . ij) behaves for small |a;| as follows: 

^ (2i)^r(o- + 1 - s) 
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with a as in (|8.2p . This function satisfies (j8.1|) and (j8.3|) . 

It remains to prove that the error term in (j8.6p does not destroy ()8.H 18. 3p . We first 
note that for r > K/\uj\ (with K as in Theorem lS.ip . the WKB estimates of Section [5] apply 
and show that (/> is well-approximated by the above Whittaker functions. Let us next show 
that for some 5 > 0, we can control the solution on the interval |a;|~^"'"'' < r < 
To this end, we introduce (similar to |1H Section 6]) the matrix 

M«,^(2iwr) VF«,^(2iwr) \ 



and the function 



Then $ satisfies the equation 



^ = A 



-1 



Again using the asymptotic formulas [H (13.1.32), (13.1.33), (13.5.1), (13.5.2)], one finds 
that I det > \io\/c, and we obtain the inequality 

\M' < p\^\ where p := , f „ . 

\uj\ r-^ 

Applying Gronwall's inequalities 

|$(r2)| < |$(ri)| exp (^j^ 

|$(r2) -$(ri)| < |$(ri)| exp / p\ \ p, 

we can easily control provided that the integral of p becomes arbitrarily small for 
small \uj\. To see this, we first note that the Whittaker functions are bounded near z = 
by c|z|5-^ and so < \z\'^-^^' (see [H (13.5.5) and (13.5.6)]). Thus 

H-l+^ " J\^\-^+S r2+2M - |a;|2M 1 + 2^' ' 

and choosing 5 < {1 + 2p)~^, the right side converges to zero as a; — > 0. 

On the remaining interval r < |c(;|^^^'^, we write the Schrodinger equation ()8.5p as 

-dr +UJ 5 V 



where we used the Heaviside function to truncate the potential in the region which is of 
no relevance here. Treating the operator on the left as the free operator, its solutions are 
given by Hankel functions (see |12l [T6] ) . A short calculation shows that the square bracket 
satisfies the condition that H?" [■ . .] is bounded uniformly in \uj\. This is precisely the 
condition which ensures the existence of the Jost solutions (see [121 Proof of Lemma 3.6]) 
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and again gives us control of the error terms. Furthermore, one sees that in the region 1 <^ 
r < lo'l"^''''', the fundamental solutions are well-approximated by the Hankel functions, 
which in turn are a limiting case of our above Whittaker solutions. This shows that the 
solution -0 of the untruncated equation (j8.5l 18. 6p has a limit as a; ^ 0, and that the 
asymptotics of the limit is the same as that of the Whittaker solutions. This justifies 
dropping the error term in (|8.6p . ■ 

Combining the last lemma with a convexity argument, we next show that the Green's 
function has a limit at w = 0. 

Lemma 8.2 For the Green's function G{u,v) 1^4-^ of the Schrddinger equation 12. 5\) . the 
limit 

lim G(u,v) 

exists and is finite. 

Proof. Using ()8.1l I8.3|l . the Green's function ()4.8p has a limit at uj = 0, 

lim G{n,v) = X | f^l^ff^] ^ ^ ^ 

a(0)9tJ^0 w{(t),(t)o) [ (P0[U>(P{V) ltv<U, 

provided that the Wronskian on the right side does not vanish. In order to show that this 
Wronskian is indeed non-zero, let us assume on the contrary that ^ is a multiple of <j)o. 
Note that for uj = 0, the potential V in (12. 6p is real and positive. Hence we can repeat 
the convexity argument in the proof of Lemma |4. II to get a contradiction, where now we 
use that (po tends to zero at infinity according to (|8.3p . ■ 

We next consider the case ojq 0. We first show that (j) has a well-defined limit 
as w — > Wo- 

Lemma 8.3 The following limit exists for any real ujq ^ and every u £M, 

lim </)(u) = (/)o('u) . 

The limiting function cj)Q is again a solution of the Schrddinger equation \2. 5\) with the 
asymptotics 113.15]) . 

Proof. We cannot introduce </>o directly via the iteration scheme (j5.5p because for real lu 
the factor exp(— 2 VV) is for large x increasing polynomially like (x — n)^''. In order 
to bypass this problem, we introduce a convergence generating factor; namely, we set 
for UJ = ujQ 

Eio) = 1 

E(^+^){u) = lim r e-'-^^^^^ll-e-'f^^\E'^'\x)dx, 
and define (po by 

oo 

^q{u) = ^S(')(n) a{u) . (8.8) 
1=1 
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Let us verify that this iteration scheme is well-defined and defines a solution of the 
Schrodinger equation (j2.5p . To this end, in (j8.7p we substitute the identity 



\-2\/vdxJ 

where, as in the proof of Theorem 15. 11 we choose p = 0,1,3 depending on whether s = |, 

1 or 2, respectively. After integrating by parts p times, the resulting integrands are 
dominated by c/x^, and thus we can take the limit 6 \ using Lebesgue's dominated 
convergence theorem. The respective estimates ()5.6l 15. 7p and <\5.9\ I5.10p for s = 1 or s = 

2 clearly remain valid for this modified iteration scheme, showing that the series (18. Sp 
converges absolutely, uniformly for sufficiently large u. 

In order to compute the u derivative of E^^^^\ we integrate by parts, take the limit 6 \ 
0, and can then compute the derivative. After this, we can re-insert the convergence 
generating factor and re-integrate by parts. This shows that in the formula for we 
may interchange differentiation with taking the limit 5 \ 0. Exactly as above, one verifies 
that the series X]i^o(-^^''')' converges again absolutely, uniformly for sufficiently large u. 
Hence (jS.Sp may be differentiated termwise, thereby showing that (/)o is indeed a solution 
of (f231) . 

Finally, to show continuity as Cs{(jJo) 3 w — > coq, we first note that because of the 
continuous dependence of the solutions of ODEs on both initial data and parameters on 
compact sets, it suffices to show continuity of <p{u) for u > ui for any sufficiently large ui. 
Again using the above integration-by-parts method, one sees that, each of the E^^\u) 
is continuous as Ce(a;o) 3 oj ^ ojq- Since for sufficiently large ui, the series converges 
absolutely, uniformly in w G Ce(a;o) U {0}, we can take the termwise limit to — > wq. ■ 

From this lemma it will follow immediately that the Green's function converges, 

hm G{u,v) = — — X <^ \ \ u [ T (8-9) 

C^(uJo)3uj^uJo w{(j),(f)o) [ 90[U)<P[V) ltv<U, 

once we have shown that the Wronskian w{(j),cj)o) is non-zero at loq. This is done in the 
next lemma. 

Lemma 8.4 For any ujq ^ 0, the Wronskian w{(j),<j)Q) ^ 0. 

Proof. Assume that w{(j),cj)o) = 0. We choose a function r] G C^([— 1, 1]). For any e < 
iOo/2, we set 

^ UJ — UJq' 



r]e{uj) = T] 



e 



and introduce the function 



R{t,u) = j dio e-''^' 7]eiu;) (l)Uu) , (8-10) 

where (j)uj{u) is the Jost solution ^. These Jost solutions have the following asymptotics, 

(j^iuiu) ~ e*'^'" as u ^ — oo 

Mu) ~ ciHu^e-^"^" + C2(cj)u-"e*'^" as n ^ oo , 
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where ci and C2 depend smoothly on u and C2(tiJo) = 0. Thus for any 5 > we can choose e 
such that 

|c2(w)| < 5 Vcj G B^{u!q) . 

Differentiating (|8.10p and using that are solutions of the Schrodinger equation (|2.5p , 
one easily verifies that R{t, u) is a solution of the Teukolsky equation (j2.ip (with the 
angular dependence separated out). Near the event horizon, R(t,u) has the following 
asymptotics, 

su su 

R{t,u) ~ ^ / dwe-*'^*%('^)e*^" = ^fis{t-u), (8.11) 
where % denotes the Fourier transform of ?7e. Similarly, near infinity, 

R{t,u) ~ —{c{7fe){t + u) + —{c^e){t-u). (8.12) 

Being the Fourier transform of a smooth function supported in ^^(ll'o), the functions fj, 
cirfi; and C2rfe all have rapid decay on the scale e~^, i.e. 

SUp|x|" (|?7e| + |cl%| + |C^|) (x) < . 



Furthermore, the function {o2rfe) is pointwise small, 



\{o^e){t-u)\ < 



similarly, all its derivatives are pointwise small. The formulas (jS.lip and (j8.12p are valid 
near u = — oo and u = oo, respectively. Since u; is in a compact set, the error terms in the 
asymptotics are bounded uniformly in time. 

The asymptotics (18.111 18.12P contradict the conservation of physical energy. Namely, 
for large negative times, (j8.1ip describes a wave of positive energy coming from the event 
horizon. However, for large positive times, the contribution of (jS.lip as well as the first 
summand in (I8.12p decay rapidly in time, whereas the energy of the second summand 
in (j8.12p . which describes a wave moving to infinity, can be made arbitrarily small by 
choosing 5 small. ■ 

We remark that if the above energy argument is made more quantitative, it even yields 
that 

ci 



C2 



< 1 . 



This is in complete agreement with the numerical result Z < 1 in the case a/M = 
obtained in [191 p. 658ff], keeping in mind that, after a time reflection, the quantities 
and Zout as introduced in p9] are multiples of our coefficients C2 and ci, respectively. 

Using Lemmas 18.21 and 18.41 for every ujq G M, we can introduce the function ^H^^q as 
the limit of the integral kernel of the resolvent from the lower half plane; namely, 

^uio{u,v) := lim R^{u,v) . 

In Theorem 17.31 we can take the limit -R — > oo and deform the lower contour £i onto the 
real axis, the upper contour ^2 onto the line Imw = gfy, to obtain the following integral 
representation, valid for all t G M. 
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Theorem 8.5 For any spin s S {^,1,2}, the solution of the Cauchy problem for the 
separated Teukolsky equation 112. 3\) with initial data = {4',dt4>)\t=o G Cq"(M)^ can be 
written as 

2m J V io + i J 

+ e-^URM{u) + ^)d., 

where both integrals are -convergent. 

Proof. As shown in the proof of Theorem 17.31 inserting the term ^o{u)/{uj + i) into the 
integrand in (j7.5p does not change the value of the hmit. According to (j7.3p . the resulting 
integrand is bounded near infinity by C/\uj\'^, and hence we can take the limit i? ^ cxd in 
the Lebesgue sense. ■ 



9 Proof of Decay 

We now prove our main theorem. 

Proof of Theorem \l.l[ As discussed in the introduction, the conservation of energy implies 
that the solution <I> of the Cauchy problem (jl.lU1.2( ) is bounded in Lf^^, uniformly in time. 
Differentiating the Teukolsky equation with respect to t, one sees that the derivatives 9"<I> 
are also solutions and are thus also bounded in Lf^^. Since the spatial part of the Teukolsky 
equation (jl.ip is uniformly elliptic away from the event horizon, we conclude that all 

9 2 2 

spatial derivatives of ^ are bounded in Lf^^. Using the Sobolev embedding H^^^ ^ ^k^c^ 
we conclude that ^ can be bounded in L"^^, uniformly in time, by a Sobolev norm of the 
initial data, i.e. for any compact set K C {ri, oo) x S^, 

sup|$(t)| < c 11(^0, ^i)||h2,2 , (9.1) 

K 

where c depends only on K and the support of the initial data. 

Decomposing the initial data into spin- weighted spherical harmonics [14], 

oo I 
l=s m=—l 

the Sobolev norm decomposes into a sum over the angular momentum modes, 

2 

^2,2 



(^0,^1)11^2 = E U^-K'^^J'!'™) 



Lm 



Since the series converges absolutely, for any e > there is an integer Iq such that 



(>Zo m=—l 



< £. 
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Hence in view of (|9.ip . the contribution of the large angular momentum modes to can 
be made pointwise small, uniformly in time. 

For the remaining finite number of angular momentum modes, we use the integral rep- 
resentation of Theorem l8.51 The integral on the real line tends to zero as f — > — oo by virtue 
of the Riemann-Lebesgue lemma. The integral on the line Imw = gfj can be bounded by 
a constant times exp(^^) and thus tends to zero exponentially fast as t — > — oo. ■ 



10 General Remarks 

We first discuss the case s = ^ of the massless Dirac equation. At first sight, it might seem 
paradoxical that in this case, the solution of the Cauchy problem has two different integral 
representations, one being the representation obtained in [9] where the w-integral runs over 
the real axis, the other being that given in Theorem l8.51 where uj is integrated over two lines 
in the complex plane. This, however, is no contradiction, as one can understand as follows. 
In [U] the massless Dirac equation is considered as a first-order system. The Teukolsky 
equation, on the other hand, is a second order equation for a single component of the 
Dirac system. It is obtained from the Dirac equation by multiplying with a particular first- 
order operator. The transformation from the Dirac equation to the Teukolsky equation 
completely changes the spectrum of the involved operators. Whereas the Hamiltonian of 
the Dirac system is self-adjoint and thus has a purely real spectrum, the Hamiltonian of the 
Teukolsky system has non-real essential spectrum (see Proposition 14. 4p . The differences 
in the integral representations reflect these differences in the spectra of the corresponding 
Hamiltonians. 

It is important to note that the differences in the spectral representations do not imply 
different long-time dynamics. To explain this better, let us consider the integral repre- 
sentation of Theorem 18.51 in the limit t — > -|-oo. In this limit the exponential factor e~*'^* 
in the second integral grows exponentially in time, suggesting that ^{t,u) should also 
increase in time. However, this reasoning is not valid because, as explained in the in- 
troduction, by considering the time-reflected Teukolsky equation for — s and using the 
Teukolsky-Starobinsky identities, we conclude that ^{t,u) indeed also decays as t — > -|-oo. 
Another way of seeing why the naive reasoning is not valid is to consider the asymptotics 
near u = +oo. Then the the fundamental solutions 0(n) appearing in the resolvent of the 
second integral, decay exponentially as ti ^ oo. This leads to an exponential damping 
of an outgoing wave moving towards infinity, which just compensates the exponential in- 
crease of the factor e"*"^*. This argument illustrates how Theorem 18 . 51 describes the correct 
dynamics in the asymptotic limit of wave packets near spatial infinity. 

We end the paper by discussing to which extent our arguments carry over to the Kerr 
metric. The first complication in Kerr is that the separation constant A depends now on lo, 
and thus the sum of all angular momentum modes must be carried along at each step, 
and this would require additional estimates to control the infinite sum. Apart from this 
additional complication, our arguments in Sections [2H7] continue to hold. In Section [8] the 
considerations before Lemma 18.41 could also be extended, provided that the sum over the 
angular momentum modes can be controlled. However, the energy argument of Lemma [8. 41 
no longer works in the Kerr geometry due to the presence of the ergosphere, where the 
physical energy density need not be positive. The numerics carried out by Press and 
Teukolsky |19j indicate that even in the Kerr geometry, the Wronskian w{(j),(j)Q) has no 
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zeros on the real line. A possible strategy to make this rigorous would be to replace our 
above energy argument by a causality argument in the spirit of ^12^ Section 7]. However, 
this would make it necessary to analytically extend the resolvent across the real line. In 
principle, this could be achieved by estimating the higher ^-derivatives of (j) similar to |12[ 
Lemma 3.4]. However, this approach seems to be technically very demanding. 

A The Energy Density of Gravitational Waves 

The Bel- Robinson tensor Qijki is defined in terms of the Riemann tensor Rijki by (see }17]) 

^ijki — -ft n^j I -\- n ii^j^ n^j i , 
where *R is the dual of the Riemann tensor given by 

*R'\, := i e'^'^'Ram , 

and e is the totally anti-symmetric tensor. Here we are concerned with perturbations of 
the Schwarzschild metric. We denote the perturbations of the Riemann tensor by AR, 

^Rijkl = Rijkl — {R^)ijkl ; 

where R^ denotes the Riemann tensor in the Schwarzschild geometry. In the Teukolsky 
framework, Ai? describes linear perturbations which are non-spherical, meaning that they 
are orthogonal to R^ in L'^{S'^). Hence Eg^ as defined by (jl.4p can be rewritten as 

^gv = / Ur'-^o ^Rao% + *^R''\o *Ai2,o'o) dfi . (A.l) 

J i=const 

This expression is quadratic in AR. Furthermore, we see from (|1.4p that Eg^ differs 
from E, p.3p . only by the constant energy of the Schwarzschild metric, and therefore Eg^ 
is conserved due to the fact that the Bel-Robinson tensor is divergence- free [17J. 

It remains to show positivity of the integrand in (|A.ip . Using that the metric is diagonal 
and that the Riemann tensor is anti-symmetric in its first two and last two indices, we 
have 

3 

AR^\o ARjo = 5°° Yl 5""/^Ai?„o/30 ^RaOfSO , 

which is obviously non-negative because g^^ > and g°^'^,g^^ < 0. Since the dual of the 
Riemann tensor has the same symmetry properties, the last summand in (|A.ip is likewise 
non-negative. 

Acknowledgments: We would like to thank Bernard Whiting for his very careful reading 
of this paper and for several valuable comments. We are grateful to the Alexander-von- 
Humboldt foundation for its generous support. We would like to thank Thierry Daude, 
Niky Kamran and Shing-Tung Yau for helpful disussions. 



30 



References 



[1] M. Abramowitz, LA. Stegun (eds) "Handbook of Mathematical Functions with For- 
mulas, Graphs and Mathematical Tables, 9th printing" Dover, New York (1972), 
www .math.sfu.ca/~cbm/ aands / sub j . htm 

[2] V. de Alfaro, T. Regge, "Potential Scattering," North-Holland Publishing Company, 
Amsterdam (1965) 

[3] A. Bachelot, "Superradiance and scattering of the charged Klein-Gordon field by a 
step-like electrostatic potential," J. Math. Pures Appl. 83 (2004) 1179-1239 

[4] P. Blue, A. Soffer, "The wave equation on the Schwarzschild metric II: Local decay 
for the spin 2 Regge Wheeler equation," ,gr-qc/031006 6, J. Math. Phys. 46 (2005) 
012502; Erratum gr-qc/0608073 

[5] J.M. Bardeen, W.H. Press, "Radiation fields in the Schwarzschild background," ,J. 
Math. Phys. 14 (1973) 7-19 

[6] S. Chandrasekhar,"The Mathematical Theory of Black Holes," Oxford University 
Press (1983) 

[7] J.M. Cohen, L.S. Kegeles, "Electromagnetic fields in curved spaces: a constructive 
procedure," Phys. Rev. D (3) 10 (1974) 1070-1084 

[8] M. Dafermos, I. Rodnianski, "A proof of Price's law for the collapse of a self- 
gravitating scalar field," |gr-qc/0309115l Invent. Math. 162 (2005) 381-457 

[9] F. Finster, N. Kamran, J. Smoller, S.-T. Yau, "The long-time dynamics of Dirac 
particles in the Kerr-Newman black hole geometry," Adv. Theor. Math. Phys. 7 (2003) 
25-52 

[10] F. Finster, N. Kamran, J. Smoller, S.-T. Yau, "Decay rates and probability estimates 
for massive Dirac particles in the Kerr-Newman black hole geometry," |gr-qc/0107094 
Commun. Math. Phys. 230 (2002) 201-244 

[11] F. Finster, N. Kamran, J. Smoller, S.-T. Yau, "An integral spectral representation of 
the propagator for the wave equation in the Kerr geometry," gr-qc/0310024^ Commun. 
Math. Phys. 260 (2005) 257-298 

[12] F. Finster, N. Kamran, J. Smoller, S.-T. Yau, "Decay of solutions of the wave equation 
in the Kerr geometry," gr-qc/0504047, Commun. Math. Phys. 264 (2006) 465-503 

[13] J. Friedman, M. Morris, "Schwarzschild perturbations die in time," J. Math. Phys. 41 
(2000) 7529-7534 

[14] J.N. Goldberg et al, "Spin-s spherical harmonics and 9," J. Math. Phys. 8 (1967) 
2155-2161 

[15] B. Kay, R. Wald, "Linear stability of Schwarzschild under perturbations which are 
nonvanishing on the bifurcation 2-sphere," Class. Quantum Cravity 4 (1987) 893-898 

[16] J. Kronthaler, "The Cauchy problem for the wave equation in the Schwarzschild 
geometry," gr-qc/0601131, J. Math. Phys. 47 042501 (2006) 



31 



[17] S. Klainerman, F. Nicolo, "The Evolution Problem in General Relativity," Birkhauser 
Boston (2003) 

[18] CO. Lousto, B.F. Whiting, "Reconstruction of black hole metric perturbations from 
the Weyl curvature," Phys. Rev. D 66 (2002) 024026 

[19] W.H. Press, S.A. Teukolsky, "Perturbations of a rotating black hole. II. Dynamical 
stability of the Kerr metric," Astrophys. J. 185 (1973) 649 

[20] R.H. Price, "Nonspherical perturbations of relativistic gravitational collapse I, scalar 
and gravitational perturbations," Phys. Rev. D (3) 5 (1972) 2419-2438 

[21] T. Regge, J. A. Wheeler, "Stability of the Schwarzschild singularity," Phys. Rev. (2) 
108 (1957) 1063-1069 

[22] S. Teukolsky, "Perturbations of a rotating black hole. I. Fundamental equations for 
gravitational, electromagnetic, and neutrino-field perturbations," Astrophys. J. 185 
(1973) 635-647 

[23] R.M. Wald, "On perturbations of a Kerr black hole," J. Math. Phys. 14 (1973) 1453- 
1461 

[24] R.M. Wald, "Construction of solutions of gravitational, electromagnetic, or other 
perturbation equations from solutions of decoupled equations," Phys. Rev. Lett. 41 
(1978) 203-206. 

[25] B.F. Whiting, "Mode stability of the Kerr black hole," J. Math. Phys. 30 (1989) 
1301-1305 

[26] B.F. Whiting, L.R. Price, "Metric reconstruction from Weyl scalars," Class. Quantum 
Grav. 22 (2005) 589-604 



Felix Finster Joel Smoller 

NWF I - Mathematik Mathematics Department 

Universitat Regensburg The University of Michigan 

93040 Regensburg, Germany Ann Arbor, MI 48109, USA 

Felix . FinsterOmathematik . uni-r . de smollerSumich . edu 



32 



